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NUMERICAL STUDY OF WAVE PROPAGATION IN A NON-UNIFORM FLOW 


ALEX POVITSKY* 

Abstract. The propagation of acoustic waves originating from cylindrical and spherical pulses, in a 
non-uniform mean flow, and in the presence of a reflecting wall is investigated by Hardin and Pope approach 
using compact approximation of spatial derivatives. The 2-D and 3-D stagnation flows and a flow around 
a cylinder are taken as prototypes of real world flows with strong gradients of mean pressure and velocity. 
The intensity and directivity of acoustic wave patterns appear to be quite different from the benchmark 
solutions obtained in a static environment for the same geometry. The physical reasons for amplification and 
weakening of sound are discussed in terms of dynamics of wave profile and redistribution of acoustic energy 
and its potential and kinetic components. For an acoustic wave in the flow around a cylinder, the observed 
mean acoustic pressure is approximately doubled (upstream pulse position) and halved (downstream pulse 
position) in comparison with the sound propagation in static ambient conditions. 

Key words, aeroacoustics, non-uniform mean flow, higher-order compact scheme, amplification of 
sound, stagnation flow, flow around cylinder 

Subject classification. Fluid Mechanics 

1. Introduction. This study investigates numerically the influence of strong mean flow gradients on 
the directivity and strength of sound waves propagating in such a flow. 

The goal of this study is to get computational insight into physical mechanisms of angular redistribution 
of acoustic energy and acoustic pressure. The considered mean flows include potential flow around a circular 
cylinder, the 2-D and 3-D stagnation flows, and a 90° corner flow. These idealized mean flows mimic real 
world flows in areas of strong sound reflection and generation such as a leading edge of a wing or a turbine 
blade, a protruding corner of a wall cavity, a wing-fuselage intersection, an impingement area of a jet, and 
a flow behind a bluff body. 

The amplification or weakening of sound propagating in non-uniform flows has attracted considerable 
attention of researchers. However, available computational studies are mainly restricted to one-dimensional 
mean flows. 

For one-dimensional flow in a duct with variable cross-section, the acoustic pressure appears to con- 
centrate towards the duct centerline or the duct wall for upstream and downstream sound propagation, 
respectively [11], [15]. Motivated by engine noise, the authors of these studies use the wave envelope method 
(WEM) to study propagation of sound in one-dimensional and quasi one-dimensional mean flows This method 
is based on consideration of the governing acoustic variables as sums of harmonics. Additional assumptions 
are made about rapid axial oscillations and slowly varying amplitude and phase corrections caused by the 
non-uniformity of the flow. 

Atassi and Grzedzinsky [2] considered propagation of unsteady disturbances in flows around bodies 
with a stagnation point. For incompressible mean flows, the aeroacoustic problem is formulated in terms 
of an integral equation of the Fredholm type. The considered streamlined body and the corresponding 
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potential mean flow result from superposition of a uniform flow and a source. Yet, additional assumptions 
and simplifications (some of them depend on the type of potential background flow) are needed to perform 
integration. Authors [2] present the disturbance pressure along the stagnation streamline that shows increase 
of the disturbance pressure for an upstream harmonic disturbance. 

Hardin and Pope approach [7] implies consideration of a flowfield as a sum of the mean velocity and 
a disturbance field and solution of disturbance variables in time and space domain by explicit integration 
in time and use of appropriate spatial discretization. The computational methodology used in the current 
study is based on above approach. 

In order to use an analytical solution for mean flow, the flow in the current study is assumed to be 
potential, i.e., inviscid and incompressible. The speed of sound (relative to the local mean flow velocity) is 
assumed to be equal to unity. These assumptions are fair unless the Mach number exceeds 0.5. 

In this paper, the relative orientation of the direction of wave propagation, mean pressure gradient, 
and the mean flow vector are studied to show the amplification or weakening of the acoustic pressure. 
The consideration starts with two-dimensional stagnation mean flows, then the approach is expanded to 
three-dimensional stagnation mean flows, and, finally, applied to the inviscid flow around circular cylinder. 

In all considered cases, the obtained acoustic pressure is compared with that for the wave propagation 
under static ambient conditions and for the same geometry of surrounding rigid boundaries. 

Variation of directivity of acoustic pressure implies changes in directivity of acoustic energy in comparison 
with the sound propagation under static ambient conditions. This paper studies angular redistribution of 
acoustic energy and its potential and kinetic components in non-uniform mean flows. The acoustic energy 
is a sum of potential and kinetic components where the potential energy is closely related to the averaged- 
in-time root mean square of acoustic pressure which is usually measured by experimentalists as an indicator 
of noise intensity. 

Our recent study [12] shows that the generation of sound waves originated from entropy sources is caused 
by baroclinic generation of disturbance vorticity as a result of interaction of the mean flow pressure gradient 
and the wave density gradient. In turn, this generated vorticity causes the generation (or sink) of acoustic 
energy by the source term in the acoustic energy equation. Order of magnitude analysis suggests that the 
source term can not be neglected in the acoustic energy equation. This study presents vorticity and energy 
source terms obtained from numerical simulation. For the propagation of acoustic waves originated from 
an acoustic pulse (as opposed to the sound generation from entropy or vorticity pulses [12]), The major 
effect of the non-uniformity of the mean flow on the wave propagation (as opposed to the wave generation 
from entropy or vorticity pulses [12]) is the angular redistribution of acoustic energy and kinetic-to-potential 
transform whereas the newly generated acoustic energy appears to be minor. 

The article is built as follows. In Section 2, the wave propagation originated from a cylindrical 2-D 
source in 2-D stagnation-type flows is studied and numerical results about amplification of acoustic pres- 
sure, alternation of wave pattern, and acoustic energy redistribution are presented. In Section 3, the wave 
propagation from a spherical source is considered in 3-D stagnation flows. The symmetric 3-D stagnation 
flow with a single stagnation point and the 3-D stagnation flow with stagnation line flow are taken as typical 
representatives of three-dimensional stagnation mean flows. In Section 4, the inviscid flow around a circular 
cylinder is considered for various subsonic Mach numbers of the mean flow. Amplification of sound (for 
upstream wave propagation) and weakening of sound (for downstream wave propagation) is studied in this 
Section. 
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2. Sound propagation in a 2-D stagnation-type mean flow. The common type of mean flow is 
that the flux across a given surface is equal to zero, either because the surface is a symmetry surface or 
because the surface is the boundary of a rigid body [3]. If two straight zero-flux boundaries intersecting at 
an angle it /n the stream-function field for incompressible flow is given by 

xp = Ar n sin(n9), (2.1) 

where (r, 8) are polar coordinates. Note that xp = 0, if 9 — 0 or 8 = 7r/n. For example, the 2-D stagnation 
flow 


U = x,V = -y (2.2) 

is a union of two flows with n = 2 where the dividing streamline is the symmetry axis of the flow. 

Introducing a disturbance, instantaneous velocities and density are considered as sums of the known 
steady incompressible mean flow (see above) and the unsteady compressible isentropic disturbance 

u = U + ut,v = V + v/,p = P + pt, p = 1 + pi. (2.3) 


The normalized speed of sound is assumed to be equal to unity and, therefore, pressure and density distur- 
bances are equal for acoustic pulse. Substituting above sums to the Euler equations, the dynamic equations 
for unsteady (disturbance) components of mass and momentum fluxes are obtained 


Qt + / * + 9y — 0 > 


(2.4) 


where 



( p! ^ 


/ (1 + pt)(U + ut) 

) 


( (1 + p/)(V + vt) \ 

q= 

(1 + p/)(U + ut) 

, / = 

(1 + pt)(U + ut) 2 + {P + pr) 


i 9 = 

(1 + pt)(U + ut)(V + vt) 


{ (l + p ,)(y + v>) ) 


\ (1 + pt)(U + ut)(V + vt) 


\ (1 + pt)(V + v/) 2 + (P + pt) ) 


(2.5) 


The equations are 
domain 


linearized and solved with respect to unknown disturbance variables in space-time 


du 

dt 

dv 

dt 

dp 

dt 


dp du dU 

dx dx dx 


,.du dU 
dy dy V 
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du dv dp dp 

dx dy dx dy 


(2.6) 


The spatial derivatives of the mean flow variables for arbitrary n are presented in Appendix A. The non-linear 
set of equations for the 2-D stagnation mean flow is presented in our study [12]. 


2.1. Numerical algorithm and the code validation. The solution is advanced in time in five sub- 
stages per time step using a low-storage explicit-in-time fourth-order Runge-Kutta (RK) scheme scheme 
proposed by Williamson [16] and implemented by Wilson et al. [17]. 

The spatial derivatives of unknown disturbance variables are approximated using compact finite differ- 
ence schemes [10] applied to natural variables ut, v> and pt 


a L J i-i + U% + a U%+\ — ^r: 


(2.7) 
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where U i , and L r (+ , are unknown derivatives at grid nodes * — l,i, and i + 1, Ax is the grid step and 

Ui- 2 ,~-,Ui + 2 denote values of variable U at grid nodes i — 2 + 2. The values of U are taken from the 
previous stage of the RK scheme. Here the classical fourth-order Pade compact scheme (a = 0.25, a = 1.5 
and b = 0) which has a tridiagonal matrix for the right and left sides of (2.7) is used. 

At all boundaries but a rigid plate characteristic inflow or outflow boundary conditions are applied 
to disturbance variables. At a rigid surface, reflection boundary conditions are used for pressure distur- 
bance and for the parallel to the surface component of velocity disturbance whereas the normal to the rigid 
surface component of velocity is taken equal to zero. Discretization of spatial derivatives in the direction 
perpendicular to a boundary are computed by one-sided finite differences at all boundaries [4] . 

The system (2.7) with near-boundary one-sided differences is solved by version of Gaussian elimination 
for narrow banded systems known as the Thomas algorithm [9] applied to each line of numerical grid in all 
spatial directions. 

To validate the code against the known analytical solution [8], the static ambient conditions (i.e., the 
mean flow velocity is equal to zero) are taken. In this case, the right-hand side of system (2.6) is simplified 
as follows 

du dp 

dt di- 

do _ dp 
dt dy 

dp du dv 

s = % < 2 - S) 

The initial 2-D cylindrical acoustic pulse is given by 


p = p = e exp 


g ( x ~ x c? + (y - yc? 

a 


(2.9) 


where e = 0.01, a = ln(2)/9,d = 60 is the normalization coefficient, (x c ,y c ) is the pulse center coordinates. 

To check the symmetry of solution, the computational domain is taken as a union of [—1,1] x [—1,0] 
and [—1, 1] x [0, 1] subdomains separated by a rigid flat plate located at y = 0. The couple of initial pulses 
(2.9) is located at (0,±0.25). Isolines of acoustic pressure for this problem at 160 x 160 grid per subdomain 
are presented in Fig. 1 at t = 0.5. The isobar values are uniformly distributed between the maximum and 
minimum of acoustic pressure. Although the Thomas algorithm is, strictly speaking, non-symmetric (as any 
Gaussian elimination), the acoustic pressure field in Fig. la is perfectly symmetric. 

The results of grid refinement study are presented in Fig. lb for [—1,1] x [0, 1] computational domain 
covered with 80 x 80, 120 x 120, and 160 x 160 numerical grids. The visible wiggles exist only for the coarse 
grid and the solution approaches the analytical solution including most critical points of the maximum and 
minimum of the reflected wave. For the transmitted wave, the numerical solution practically coincide with 
the analytical solution [8] for finer grids. Therefore, the grid 160 x 160 has been adopted in this study. 

2.2. Distribution of acoustic pressure. To study the acoustic pressure distribution in a non-uniform 
flow, the 2-D stagnation inviscid mean flow is taken as an example (Case 2- A). The velocity field of the 
stagnation mean flow and isobars of the initial acoustic pulse, centered at (0,0.25), are presented in Fig. 2a 
(left). The speed of sound is equal to unity and is presented in Fig. 2a for scale. As far as the stagnation 
flowfield satisfies Eq. (2.2), the Mach number at a point is equal to the local distance from the origin, i.e., 
the Mach number at the location of the center of initial pulse is equal to 0.25. 
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The propagation of sound is described by the system (2.6) where the mean flow is given by (2.1) with 
n = 2. Grid refinement study for Case 2-A is presented in [12] and shows that the results on the 160 x 160 
and 200 x 200 grids are quite closed to each other that confirms the use of the former grid for simulations. 

The isobars of acoustic pressure are presented in Fig. 2a (right) at t = 0.5. Approximate centerline 
coordinates of the maximums of acoustic pressure for the transmitted and reflected waves are 0.5 and 0.25, 
respectively. Note that the Mach numbers at these points are equal to their y— coordinate values (see above). 
As expected, the acoustic isobars obtained in the static ambient conditions, denoted here as Case 2-B, (Fig. 
la) are different from those in the Case 2-A at the same time moment. In the Case 2-A, the y — coordinate 
of the wave front at the neighborhood of the centerline is smaller and the x — coordinate of the wave front 
near the wall is larger than those in the Case 2-B. The explanation is straightforward and is based on the 
fact that the speed of sound in a steady frame is equal to c -I- U, where c is the local speed of sound and U 
is the local velocity of the flowfield. 

The density of isobars is larger at the neighborhood of the centerline and smaller near the wall in the 
Case 2-A in comparison with those in the Case 2-B. In Fig. 3a, the acoustic pressure profiles at the centerline 
are presented for both cases at t = 0.5. The maximum absolute value of acoustic pressure is approximately 
80% larger for the transmitted wave and two times larger for the reflected wave than those in the Case 2-B. 
Recall, that the speed of wave propagation is c + U. While the wave propagates upstream the decelerating 
subsonic flow, its front moves slower than its back. The situation is opposite for the accelerating flow (which 
is the case near the wall). Thus, the wave front becomes narrower at the neighborhood of the centerline and 
wider while the stagnation flow spread along the wall in comparison to those in static ambient conditions. 

Accelerating mean flow (backward 2-D stagnation flow) is considered to study the influence of the 
mean flow direction on the wave patterns (Fig. 2b). The initial position of acoustic pulse, the mean flow 
streamlines and the mean pressure distribution remain the same as for the Case 2-A. (compare Fig. 2b with 
Fig. 2a). However, the direction of the maximum acoustic pressure has been considerably changed from 
being perpendicular to the wall to being parallel to it. Still, the maximum acoustic pressure corresponds to 
the direction where the sound propagates upstream of the mean flow. 

To study the influence of the background flow while the acoustic pulse is located off the centerline and 
two reflecting surfaces are present, the 90°— corner geometry is considered (Fig. 7). The rigid plains coincide 
with the Cartesian axis. The potential flow in this geometry is described by a single stagnation flow with 
n — 2 flowing to the right. Initial position of acoustic pulse is taken (0.25,0.25), where the local Mach 
number of the mean flow is approximately equal to 0.35. 

Computational results in terms of acoustic pressure are compared with those obtained in the Case 2-B. 
In the latter case the maximum acoustic pressure is in the bisector direction and the cylindrical shape of 
wave pattern is symmetric with respect to the bisector (Fig. 7b). In presence of the mean flow, the waves 
flatten, especially in the neighborhood of the y— axis, where the sound propagates upstream of the mean 
flow (Fig. 7a). The maximum acoustic pressure increases about 40% in comparison with that in the static 
ambient conditions. 

2.3. Distribution of acoustic energy. To get physical insight into amplification of acoustic pres- 
sure, acoustic energy is considered here. Transport equation of acoustic energy density for irrotational and 
incompressible mean flow and adiabatic acoustic disturbance is given by [6]: 

dE 

— — b V/ = uf ■ (U x c or), (2.10) 
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where 


E = p/ 2 / 2 + uf 2 /2 + p/uf ■ U 


is the acoustic energy density and 


I = (pi + uf ■ U)(uf + p/U ) 


( 2 . 11 ) 


(2.12) 


is the acoustic energy flux. 

The acoustic energy is a sum of potential energy Pot = p/ 2 / 2, that mimics the strength of acoustic 
pressure in compressible media, and kinetic energy K = u/ 2 / 2 + p/ui ■ U. 

The centerline distribution of E is presented in Fig. 3b, where the curve 1 corresponds to the Case 
2- A and the curve 2 corresponds to the Case 2-B. For the reflected wave, the local maximum of E in the 
Case 2-A is about three times greater than that in the Case 2-B. Its location is the same as that for the 
maximum acoustic pressure (compare Fig. 3b to Fig. 3a) and the ratio of acoustic energies is similar to the 
ratio of squares of acoustic pressures (four times). It is not the case for the transmitted wave, where the 
maximum values of E are almost the same in cases 2-A and 2-B in spite of the obvious difference in the 
acoustic pressure values. 

To show the total acoustic energy along the centerline, the integrals 

Ie= r E(y)dy (2.13) 

Jo 


and 


Ip = [ Pot(y)dy 
Jo 


(2.14) 


are plotted as functions of y at the time moment t = 0.5 in Fig. 3c. Obviously, the integrals reach their 
maximum when the coordinate y exceeds the front coordinate of the transmitted wave. The values of these 
maximums for integral of acoustic energy differ only by 15% whereas the difference in integrals of potential 
energy is about two times. Thus, the amplification of sound along the centerline occurs mainly by increase 
of the potential part in total acoustic energy. 

On the contrary, for the bisector direction (45°) and for a near-wall direction (75°) the part of potential 
energy remains approximately the same for the stagnation mean flow and for the static ambient conditions. 
Therefore, the weakening of sound near the wall occurs mainly due to the angular redistribution of acoustic 
energy. The direction-dependent behavior of acoustic energy is explained by the relative value of the last 
term in (2.11). If sound propagates upstream of the mean flow, the vectors in the inner product ui ■ U are 
collinear, therefore, the variation of U in a non-uniform mean flow leads to more prominent reduction of the 
kinetic energy in comparison with other directions. The kinetic energy transforms into potential energy and 
leads to increase of potential part in the acoustic energy. 

To study the generation of acoustic energy by means of interaction of the wave with the non-uniform 
mean flow, we present here the simplified form of the vorticity transport equation [12]: 


— = Vp/ x VP. (2.15) 

Note that the presence of the pressure gradient in the mean flow (i.e., non-uniformity of the flow) leads 
to non-zero right-hand side of above equation. Otherwise, the vorticity wi remains equal to zero as in an 
irrotational mean flow. In turn, the non-zero c of leads to the non-zero source term in the energy equation 
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(2.10). The vorticity is generated where the acoustic wave passes, i.e. Vpf ^ 0. As far as ut is a multiplier 
in the right-hand side of Eq. (2.10), the exchange of energy between the steady mean flow and the unsteady 
disturbance occurs only in presence of the wave. The generated vorticity flows with the local speed of 
background flow and eventually the vorticity is left behind the wave, that moves with the sonic speed. The 
instantaneous pattern of the vorticity generation source and the acoustic energy source at t = 0.5 are shown 
in Fig. 5. As expected, areas of maximum acoustic energy generation coincide with instantaneous positions 
of the reflected and transmitted waves. 

Approximately, the maximum value of the energy source term is located at the bisector direction (45°). 
In Fig. 6, the acoustic energy and its source are shown as functions of the distance from the origin along the 
bisector line at t = 0.5. The maximum value of the source is approximately one order of magnitude smaller 
than that for the acoustic energy. Therefore, the local gain of acoustic energy is small in comparison with 
the angular redistribution of acoustic energy and kinetic-to-potential transforms considered above. 

3. Acoustic pulse in the 3-D stagnation flow. In this section we consider a spherical pulse prop- 
agating in a 3-D stagnation flow. The 3-D stagnation flow with a single stagnation point at the origin is 
denoted as Case 3-A and given by 


W = —z, U = 0.5a:, V = 0.5 y. 
The stagnation flow with the stagnation line (Case 3-B) is 

W = -z, U = x, V = 0. 


(3.1) 


(3-2) 


The developed code is readily available for 3-D parallel computations as well [13]. 

In Case 3-A, the propagation of disturbance is described by the following system of equations 


du 

dt 

dv 

dt 

dw 

Hi 

dp 

dt 


— "> + *& 

dp „ dw dw , dw 

-d-z + °- 5 ( - y d^- X d^ ) + W - Zp + Z di’ 


du 

dx 


dv 

dy 


£ +0 - 5 <- 


dx 


dp . dp 


In Case 3-B, the governing system of equations is given by 


du 

dt 

dv 

dt 

dw 

~dt 

dp 

dt 


dp 


du 


du 


x- u - xp+ z— 

dx dx dz 

dp d v _ d v 

dy X dx ^ ^ dz’ 


dp 


dw 


dw 


~a + w - zp + z—, 

dz dx dz 

du dv dw dp _ dp 

dx dy dz X dx+ ° dz 


(3.3) 


(3.4) 


The initial 3-D spherical acoustic pulse is given by an expression similar to (2.9). 

To validate the code, the computations are performed on a set of numerical grids where the wave spreads 
under static ambient conditions and in presence of the rigid reflecting plate (Case 3-C). The analytic solution 
for acoustic pressure in case of the infinite domain is given by 

P = ^(( r - *) ex P[ -a ( r - ^ + ( r + *) exp[-a(r + f) 2 ]}, (3.5) 


7 



where r is the location of the pulse center at t = 0. Solution for the semi-infinite domain is obtained by use 
of the image pulse located at — r : 


p = p(r) +p(-r). (3.6) 

Results of computations on 80 3 , 120 3 , and 160 3 grids are compared to the above analytic solution. 
Profiles of the centerline acoustic pressure on the above set of grids at t = 0.5 are presented in Fig. 8. The 
observations from the grid refinement study are similar to those obtained for the 2-D case (Section 2). The 
numerical solution on the 160 3 numerical grid practically coincides with the analytic solution and this grid 
is used for numerical simulations. 

To compare acoustic fields in Cases 3-A, 3-B, and 3-C, the isolines of acoustic pressure at t — 0.5 are 
shown in Fig. 9. In Cases 3-A and 3-C (Fig. 9a, b), the acoustic field is uniform with respect to the direction 
in the x — y plane. The presence of the stagnation flow leads to amplification of sound at the neighborhood 
of the centerline. In Case 3-B, the wave pattern is direction-dependent: the maximum elongation of the 
wave pattern along the rigid plane occurs in the section x — z, whereas the minimum elongation is in the 
section y — z. In the latter case, the elongation is closed to that for the Case 3-C (compare Fig. 9a to Fig. 
9d). Yet, the centerline position of the pulse is the same in Cases 3-A and 3-B (see Fig. 9b, c, and d). 

The centerline acoustic pressure for Cases 3-A and 3-B is presented together with its counterpart for the 
Case 3-C at t = 0.5 in Fig. 10. In spite of the different wave pattern in Cases 3-A and 3-B, the centerline 
profiles coincide for these cases. The amplification of acoustic pressure along the centerline is controlled 
by redistribution of the potential and kinetic parts of acoustic energy in this direction as opposed to the 
angular redistribution of acoustic energy (see the previous section). As far as the centerline profiles of the 
mean velocity and pressure are the same in Cases 3-A and 3-B, the profiles of acoustic pressure appear to 
be the same along the centerline. Note that this is quite a special situation where both the initial conditions 
of pulse and the background flow velocity are the same along the streamline for two different flows. 

The maximum acoustic pressure for the transmitted wave is two times larger in Cases 3-A and 3-B than 
that in the Case 3-C. The maximum acoustic pressure for the reflected wave becomes about three times 
larger in presence of stagnation flow than that for the static ambient conditions. 

The amplification of sound in the 3-D case is larger than that in the 2-D case. Recall that in both cases 
the background flow velocity and pressure are the same along the centerline. To get qualitative explanation 
of such a difference, the analytic solutions for acoustic pressure at the centerline (static ambient conditions) 
is presented for the 2-D and 3-D cases in Fig. 11a. The spatial derivative of acoustic pressure along the 
centerline is shown in Fig. lib. By inspection of Fig. 11a. and by comparison of derivatives in Fig. lib, one 
can conclude that the back and front wave fronts are sharper in the 3-D case than those in the 2-D case. 
The front of the wave propagating upstream move slower than the back of it, that leads to the amplification 
of sound (see the previous section). The sharper the wave front is, the more prominent is the phenomenon. 
Therefore, the 3-D pulse in a 3-D stagnation flow amplifies larger than that in the 2-D case. 

4. Aeroacoustics of the flow over a circular cylinder. The propagation of sound waves from 
a single acoustic pulse in the 2-D, inviscid irrotational, incompressible mean flow over an infinite circular 
cylinder is considered here. The computational domain is shown in Fig. 12. The circle of unit diameter 
(. R cy i = 0.5) is the inner boundary and a circle of 3.5 diameter is the outer boundary. At the outer boundary 
the mean flow is practically uniform. The streamlines of the mean flow and the isobars of mean pressure are 
computed analytically for potential flow [1] and presented in Fig. 12. The pulse is located at the centerline, 
in some distance from the rigid cylinder. While the mean flow is directed from left, to right, the acoustic 
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waves, originated from the pulse, propagate upstream of the mean flow. Otherwise, the acoustic waves 
propagate downstream of the mean flow. The mean flow in front of the cylinder can be described by the 
incompressible potential flow model fairly well up to M = 0.5. The after-body flow can be described by this 
model when M < 0.4. It has been shown that the flow streamlines around cylinder are quite similar to those 
for incompressible flow if M < 0.4. [14]. For higher Mach number, strong recirculating zones appear at the 
rear part of the cylinder and the shock waves take place at the mid-section of the cylinder. 

In polar (r, 9) coordinates, the equations of propagation of small disturbances in a known potential flow 
are given by 


du 

dt 

dv 

at 

dp 

dt 


dp TT du dU TT dU 1 /Tr du dU Tr dU , 1 /T , 2 

dr 1 dr dr U 1 dr P r V d6 + d6 V + V 89 p) r p(F Vv ^ 

~ u P~ir u ~ u ir p ~ + V 'S p) _ ~( puv + Uv + Vu) ' 

r dO r d6 dr dr dr rdO dO r 

_du _ dp _u_pU_ 1 dv dp 
dr dr r r r K d$ d9 h 


(4.1) 


where U, V, u, and v are radial and angular components of the mean flow velocity and disturbance velocity, 
respectively. The mean flow velocity and its derivatives are computed by [1]. The initial acoustic pulse is 
given by (2.9). The uniform grid 360 x 360 in polar computational space is used. 

In Fig. 13, isobars of acoustic pressure at t = 0.5 are presented for upstream propagation of acoustic 
pulse (Case CYL-A) at M = 0.5, downstream propagation of acoustic pulse (Case CYL-B) at M = 0.4, and 
static ambient conditions (Case CYL-C). In Fig. 14, the centerline pressure distribution is presented for 
these cases at time moments t = 0.5, 1.0, 1.5,2. In comparison with the Case CYL-C, the acoustic pressure 
is amplified at the neighborhood of the centerline in the Case CYL-A and is weakened in the Case CYL-B. 

To get the sound directivity at different Mach numbers for Cases CYL-A and CYL-B, the root mean 
square of acoustic pressure is calculated by 



where p is the acoustic pressure (p is equal to p) and T is the time period of summation. This time 
period is chosen constant and equal for all numerical experiments presented below. For temporal numerical 
integration, the integral in (4.2) is computed as sum of values of squares of acoustic pressure at all time steps 
over the period, T. The p 2 as a function of time for Cases CYL-A, CYL-B, and CYL-C at the centerline 
location —5 R cy i is presented in Fig. 15. Curves have two stairs corresponding to passing transmitted and 
reflected waves. The time period, T, covers the pass of both waves. The p rms represents the time-averaged 
potential part of acoustic energy passing through a fixed point in space. 

In Fig. 16, the p rms is shown as function of angle a from the centerline while the radial coordinate is 
equal to —5 R cy i. The curves are presented for various Mach numbers from M = —0.4 to M = 0.5 with step 
0.1. The pressure pulse is located at the distance 0.25 R cy i from the cylinder (Fig. 16a, b) and at the distance 
0.5 R cy i (Fig. 16c, d) at t = 0. The positive Mach numbers correspond to the Case CYL-A (Fig. 16a, c), the 
negative Mach numbers represent the Case CYL-B (Fig. 16b, d), and the zero Mach number is the Case 
CYL-C. The latter case is presented in (Fig. 16a-d) as a benchmark for comparison. The presented results 
are normalized in such a way that the p rrns at the centerline in the Case CYL-C is equal to unity. 

The Prms increases with the Mach number and reaches 2.2 while M = 0.5 (Fig. 16a). Even far apart 
from the centerline ( a = 60°), the p rms for M = 0.5 is about twice as many as that for the Case CYL-C 
(M = 0) in the same direction. 
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For the Case CYL-B, the p rms reduces with the absolute value of the Mach number and becomes 0.65 
at the centerline while M = —0.4. The level of p rms is smaller than that for the Case CYL-C up to a = 45°. 

The degree of increase and reduction of acoustic pressure is slightly bigger when the pulse is located more 
close to wall (compare Fig. 16a, b with Fig. 16b, d). Note, that the difference between acoustic pressures 
is quite small for these pulse locations in Case CYL-C. The reason for such a difference for non-zero Mach 
numbers is that the reflected wave propagates in the direction opposite to the far-field sound before it hits 
the wall. For instance, in the Case CYL-A the reflected wave propagates downstream before its reflection 
from the cylinder and its acoustic pressure decreases. 

5. Conclusions. The behavior of acoustic waves originating from the single pulse and propagated in 
the non-uniform incompressible background flow is studied by numerical simulation using Hardin and Pope 
approach. Higher-order compact spatial finite differences and Runge-Kutta temporal integration are used. 
Results are compared to those obtained for the same geometry in the static ambient conditions. 

It is shown, that the acoustic pressure is almost doubled while the waves propagate upstream of the 
stagnation flow and the background flow reaches M = 0.5. Alternation of the direction of stagnation flow 
change the direction of sound propagation from perpendicular to a wall to parallel to it. Sound propagation 
in a corner while in the presence of the background inviscid flow leads to amplification of sound and flattening 
of acoustic waves. 

In terms of acoustic energy, the modified sound directivity in the presence of stagnation flow is mainly 
caused by redistribution of potential and kinetic components of acoustic energy (while sound propagates 
upstream) and by angular redistribution of acoustic energy (while directions of the mean flow and sound 
propagation are far from collinear). The pump of acoustic energy from the background flow by means of 
baroclinically generated vorticity is minor. 

The propagation of acoustic wave originated from spherical pulse in 3-D stagnation flows was considered 
and the wave pattern in the neighborhood of the single stagnation point and the stagnation line are discussed. 
The amplification of sound by a mean flow is more prominent in the 3-D case than in the 2-D case because 
of steeper acoustic wave bounds. 

Finally, the propagation of acoustic pulse in the flow around the 2-D circular cylinder is modeled. 
The time-averaged root mean square of acoustic pressure (rms) is presented as function of angle from the 
centerline. To compare with the static ambient conditions, the rms of acoustic pressure at the centerline 
is approximately doubled (upstream pulse position, M = 0.5) and halved (downstream pulse position, 
M = 0.4). The amplification (upstream pulse position) and weakening (downstream pulse position) of sound 
holds for a large angular sector apart from the centerline. 

This study can be expanded to cases where a compressible mean flow is computed by a CFD method and 
the acoustic source distribution is extracted from appropriate turbulence modeling. Investigation of sound 
amplification is expected to be important for higher Mach number flows where the velocity and pressure 
gradients in the mean flow are larger than those in the low-Mach-number flows considered here. 
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Appendix A. Mean-flow variables for a flow between two straight streamlines. 

In Cartesian coordinates, the stream-function is given by 


ip = A(x 2 + y 2 )°' 5n sin(narctan(x/y)), (A.l) 

The velocity field (U(x,y),V(x,y)) is computed by differentiation of the stream-function defined above 

(A. 2) 

Using Mathematica [18], the following expressions for velocities and their derivatives are obtained 


tt, \ dtp 84 > 

u( - x ’»)=ai’ v(x ’ y) = Si 


U(x,y) = Anr n 2 (xcosn8 + ys'mnO), 
V(x,y) = Anr n ~ 2 (y cos n0 — xsin nO), 


9U(x,y) 
9U(x,y ) 

dy 


= Anr n 4 ((n — l)(t/ 2 — x 2 ) cosnO + 2 (n — l)xysmn8), 
= Anr" _4 ((n — l)(a; 2 — y 2 ) sin nd — 2(n — l)xy cosnO), 


where r = \f(x 2 + y 2 ), and 8 = arctan (x/y). 

For an incompressible and irrotational mean flow, 


(A.3) 

(A.4) 

(A. 5) 
(A. 6) 


dU(x,y ) 
dx 


8V{x,y) 8U (#, y) _ 8V(x,y) 
dy ’ dy dx 


(A. 7) 
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Fig. 2. pulse propagation in the 2-D stagnation flow: a) forward stagnation flow , A = — 0.5, and b) backward 

stagnation flow, A = 0.5, Eg. (2.1). The left column shown the mean flow and the isobars of acoustic pulse at t = 0, the right 
column presents acoustic isobars. 








Fig. 3. Acoustic pressure, acoustic energy (2.10), and integral of acoustic energy (2.13 and 2.14) along the centerline. 
Here (1) is the Case 2-A and (2)-is the Case 2-B. In graph (c), solid lines denote integral of acoustic energy (2.13), and 
dashed lines denote integral of potential part of acoustic energy ( 2 . 14 ) 
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Fig. 5. Generation of acoustic energy (Case 2- A): a) generation of vorticity, and b) generation of acoustic energy 


Fig. 6. Acoustic energy and its source along the bisector a = 45° at t = 0.5, where 1) is the acoustic energy (solid line), 
and 2) is the source of energy (dashed line). 1/(45°) is the distance from the origin along the bisector 












b. 

Fig. 11. Acoustic pulse propagation and reflection in static ambient conditions (analytic solution): a) centerline acoustic 
pressure, b) the spatial derivative of acoustic pressure along the centerline. The direction of the x— coordinate coincides with 
the centerline, graphs (1) and (2) denote 2-D cylindrical pulse (Case 2-B) and 3-D spherical pulse (Case 3-C), respectively. 
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Fig. 12. Potential non-lifting flow over a circular cylinder. Here P denote initial position of acoustic pulse, and thick 
lines denote boundaries of cylindrical computational domain. Isolines of mean pressure and streamlines are shown. 
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Fig. 14. Centerline acoustic pressure for pulse near a circular cylinder: Conditions for graphs (a-c) are the same as in 
the previous figure. 
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